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ABSTRACT 

Internal  filamentary  glass  damage  caused  by  high  power  Q-svitched 
pulse  lasers  has  been  reported  in  the  literature.  Also,  filamentary  trapping 
in  liquids  has  been  observed.  This  report  presents  a  theory  of  acoustic  beam 
trapping.  In  the  theory,  electrostriction  is  the  sound  wave  driving  force  and 
sound  wave  compressions  cause  focusing  of  the  light  wave  fields.  The  theo¬ 
retical  trapping  thresholds  are  calculated  from  the  laser  wavelength  and  the 
density,  refractive  index.  Young's  modulus,  and  Poisson's  ratio  of  a  solid 
material,  or  the  density,  refractive  index,  and  speed  of  sound  of  a  liquid 
medium.  The  acoustic  trapping  thresholds  agree  with  experimental  glass 
damage  thresholds  to  within  experimental  accuracy,  and  they  vary  the  same  way 
with  initial  beam  size.  Trapping  movies  showing  tbt  formation  of  strongly 
focused  regions  are  presented.  An  explanation  is  given  for  most  of  the  salient 
features  observed  in  the  damage  phenomenon. 

An  analysis  of  Kerr  effect  trapping  is  also  presented  for  purposes 
of  comparison. 
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SECTION  I 

PREFACE 

The  goal  of  the  research  reported  here  is  to  continue  the  develop¬ 
ment  of  a  theory  of  acoustic  beom  trapping.  The  theory  is  a  possible  explana¬ 
tion  of  the  mechanism  for  internal,  filamentary  glass  damage  by  lasers.  It 
may  also  be  an  explanation  for  some  of  the  beam-trapping  phenomena  observed  in 
liquids  when  they  are  traversed  by  a  high  intensity  laser  pulse. 

We  began  developing  the  theory  in  1965.  Sections  II  and  III  pre¬ 
rent  theoretical  results  obtained  by  August  19b6.  These  results  were  origi¬ 
nally  presented  in  the  proposal  for  this  contract.  Since  that  time,  some  im¬ 
provements  have  been  incorporated  in  Sections  II  and  III.  The  trapping 

2 

threshold  was  reduced  by  a  factor  of  n  /4  from  the  formula  in  the  proposal. 

o 

Also,  new  information  from  the  manufacturer  of  our  glass  samples  has  been  ac¬ 
quired.  The  graph  comparing  the  theoretical  and  experimental  damage  thresholds 
for  three  optical  glasses  is  now  based  on  the  new,  more  accurate  values  of  the 
glass  bulk  parameters.  The  agreement  between  theory  and  experiment  is  greatly 
improved. 

Section  IV  covers  an  analytical  solution  to  the  beam-trapping  equa¬ 
tions  for  the  case  in  which  Kerr  effect  trapping  is  the  dominant  trapping  mach- 
anism  and  where  electrostrictively  driven  sound  waves  are  weak  or  absent. 

Section  V  presents  some  of  our  computer  results,  depicting  the 
acoustic  beam-trapping  phenomenon.  An  explanation  of  the  salient  features  of 
the  glass  damage  phenomenon  is  given. 

Section  VI  presents  a  brief  summary  of  our  work  and  our  plans  for 
the  second  half  of  the  contract. 
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SECTION  II 

INTRODUCTION 


2.1  A  NEW  POWER  LIMIT 

As  higher  and  "Lgher  power  lasers  are  developed,  basic  limitations 
in  the  power-transmitting  capability  of  materials  and  propagation  media  are  be¬ 
ing  discovered.  Examples  of  limitations  are  electric  breakdown  and  beam  insta¬ 
bility.  During  our  study  of  gain  saturation  and  other  anomalies  in  stimulated 
Raman  effect  and  in  our  experimental  work  in  laser  damage  to  glass,  we  have 
identified  a  new  kind  of  beam  instability,  acoustic  beam  trapping. 

Acoustic  beam  trapping  is  caused  by  the  focusing  action  of  electro- 
strictively  driven  acoustic  waves.  These  waves  alter  the  index  of  refraction 
of  the  medium  by  the  density  changes  they  cause.  All  optical  materials  experi¬ 
ence  these  electrostriction  forces  and  photoelastic  effects. 

Acoustic  beam  trapping  sets  upper  limits  to  beam  power  which  can  be 
transmitted  in  important  materials  such  as  glass  and  air.  Furthermore,  the 
scaling  laws  for  acoustic  trapping  (power  threshold  versus  beam  size  and  pulse 
duration  time)  differ  from  the  scaling  laws  for  other  types  of  beam  instability 
such  as  Kerr  effect  or  anomalous  dispersion  trapping.  In  fact,  for  many  trans¬ 
parent  materials,  there  is  a  beam  size  and  pulse  duration  for  which  the  acousti 
trapping  threshold  is  lower  than  the  thresholds  fjr  other  known  instabilities. 

2.2  THRESHOLD  PREDICTABLE 

At  the  present  time,  we  can  predict  the  acoustic  trapping  power 
threshold  in  glassy  materials  from  a  knowledge  of  material  properties,  wave¬ 
length,  pulse  duration  time,  and  beam  size.  There  is  a  critical  power  level 
for  each  material  which  can  be  calculated  from  tabulated  material  properties 
and  the  laser  wavelength. 
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2.3  RECONCILABLE  WITH  STEAD/  STATE  TRAPPING  THEORY 

We  have  also  shown  that  for  any  given  beam  site  and  type  of  material, 
there  is  a  maximum  average  rate  of  Increase  of  power  which  can  be  transmitted 
without  causing  acoustic  beam  trapping.  So  long  as  power  is  added  to  a  beam  at 
a  slower  rate  than  this,  it  is  possible  in  principle  to  reach  the  steady  state 
trapping  threshold  predicted  by  Chiao,  Garmire,  and  Townes  (Ref.  1). 


2.4  SCALING  LAWS 

Figure  1  illustrates  the  scaling  laws  for  three  different  effects 
which  limit  the  ability  of  an  optical  material  to  transmit  a  laser  beam.  In 
the  graph,  total  beam  power  is  plotted  versus  beam  radius  on  log- log  scales. 

The  three  effects  are  electric  field  breakdown,  Kerr  effect  t* applng,  and  acous¬ 
tic  trapping. 

For  electric  field  breakdown,  as  in  laser-induced  "air  sparks,  " 

2 

Power  Threshold  _  Z'  Constant  for 'N  f  Beam  ^ 
for  Breakdown  “  V.  Material  J  ^Radius/ 

On  the  log-log  plot,  Lais  power  thresnold  is  a  line  with  a  slope  of  2.  Thus, 
electric  field  breakdown  limits  the  material  to  transmission  of  power  levels 
and  beam  sizes  in  the  right-hand  portion  of  the  graph. 


For  Kerr  effect  trapping,  as  shown  by  Chiao,  Garmire,  and  Townes 

/Ref.  1) 


Power  Threshold  for 
Kerr  Effect  Trapping 


(Wavelength) 


(Constant  for  'N 
Material  J 


Since  the  threshold  is  independent  ot  beam  size,  it  appears  as  a  horizontal  line 
on  the  graph.  For  transmission  without  Kerr  effect  trapping,  the  beam  power 
level  must  be  below  the  Kerr  effect  threshold. 


Together,  Kerr  effect  trapping  and  electric  field  breakdown  limit 
beam  transmission  to  beam  sizes  and  power  levels  in  the  lower  right  portion  of 
the  graph.  Both  of  the  effects  are  virtually  independent  of  the  pulse  duration 
time. 
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Figure  1.  Scaling  Laws  for  Kerr  Effect  Trapping, 
Breakdown,  and  Acoustic  Trapping. 
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The  acoustic  trapping  threshold  can  be  calculated  as  in  Section  III 
as  long  as  the  laser  pulse  duration  time  is  less  than  the  time  required  for 
sound  to  cross  the  beam  radius.  The  scaling  luw  for  acoustic  trapping  is 

o 

Power  Threshold  for  ,  .,v?.  f  Con'.tant  for  ^  (Beam  Radius)' 

Acoustic  Trapping  =  (Wavelength)  ^  Material  )  7p 

for  the  limiting  case  of  a  very  short  pulse  or  a  very  large  beam.  7f  tha  pulse 
duration  is  held  constant,  the  threshold  power  can  be  plotted  versus  beam  size 
as  shown  in  the  graph.  The  acoustic  trapping  power  threshold  curve  flattens 
out  at  the  bottom,  at  the  critical  power  level.  This  critical  power  level, 

Pc,  is  characteristic  only  of  the  material  and  the  wavelength. 

For  the  domain  in  which  the  pulse  duration  is  longer  than  the  time 
T  required  for  sound  to  cross  the  beam,  the  acoustic  trapping  threshold  remains 
at  P£.  This  power  level  can  be  exceeded  only  if  the  power  does  not  Increase 
more  rapidly  than  Pc/T  in  any  time  T. 

If  the  pulse  duration  is  increased  by  a  numerical  factor,  and  the 
beam  radius  is  increased  by  the  same  factor,  the  power  threshold  remains  con¬ 
stant.  Thus,  in  the  leg- log  graph,  the  curve  is  shifted  to  the  righc  by  the 
log  of  the  numerical  factor  by  which  the  beam  radius  is  increased.  Thus  the 
graph  shows  three  extra  curves  shifted  by  factors  of  10,100,  and  1000. 

In  some  materials  the  Kerr  effect  trapping  threshold  is  lower  than 
the  critical  power  level  for  acoustic  trapping.  In  these  materials  Kerr  effect 
trapping  will  occur  before  acoustic  trapping.  However,  roost  common  optical 
materials  have  a  critical  power  level  lower  than  the  Kerr  effect  trapping  thresh 
old.  For  these  materials  there- will  always  be  a  domain  of  beam  size  and  pulse 
duration  time  in  which  acoustic  beam  trapping  sets  the  maximum  transmit table 
power  for  the  material > 

2.5  TYPICAL  RESULTS  FOR  GLASS 

Acoustic  trapping  is  an  important  cause  of  laser  damage  to  optical 
glass,  as  shown  in  experiments  performed  by  Steinberg,  Atwood,  Lee,  and  Ward 
(Ref.  2).  The  theoretical  trapping  threshold  is  compared  with  the  experimental 
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damage  threshold  In  Figure  2.  Generally,  the  trapping  threshold  la  below  the 
experimental  damage  threshold  points,  as  expected.  In  particular,  note  the 
experimental  damage  threshold  curve  for  dense  flint.  The  left,  lower  portion 
has  a  slope  of  2.  In  that  region,  the  cause  of  damage  Is  probably  electric 
breakdown.  The  right,  upper  portion  fits  the  curve  for  acoustic  trappln  to 
within  the  experimental  repeatability,  and  it  scales  the  same  way.  For  fused 
silica  and  BK-7  the  agreement  between  the  acoustic  crapping  threshold  and  the 
measured  damage  threshold  is  even  better. 

The  experimental  results  are  not  attributable  to  Kerr  effect  trap- 
Ping  because  of  the  dependence  on  beam  site.  Also,  In  gla.s,  Kerr  effect  is 

so  weak  that  the  power  threshold  for  Kerr  effect  trapping  is  /,  negawatts,  which 
is  above  the  top  of  the  graph. 
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SECTION  III 

DERIVATION  OF  THE  ACOUSTIC 
BEAM  TRAPPING  THRESHOLD 

3.1  THE  DRIVING  FORCE 

Electrostriction  is  the  force  exerted  by  an  electric  field  on  a 
material  medium,  when  the  force  is  proportional  to  the  square  of  the  field. 

The  net  body  force  is  th°..  proportional  to  the  power  intensity  gradient,  in 
a  lossless  medium.  The  -elevant  permittivity  ratio  at  optical^frequencies  is 
the  square  of  the  refractive  index.  Thus,  the  net  body  force  f  per  unit 
volume  due  to  the  light  beam  is  (Ref.  3) 

7  =  (:/6nQC)  (n^  +  2)  (n^  -  1)  7  I. 

Thus  a  cylindrical  beam  will  drive  a  radially  proprgating  -ound  wave. 

3.2  THE  FHOTOELASTIC  EFFECT 

The  refractive  index  change  &n  due  to  small  acoustic  compressions 
is  defined  as 

An  =  n  -  nQ  =  (p  -  po> 

Assuming  constant  polarizability  per  molecule,  may  be  calculated  by 

differentiating  the  Clausius-Mosotti  relation.  Then  we  obtain 

An  =  (1/6  nQ) (n0  +  2) (no  -  1)  O 

where  a  is  the  normalized  compression. 

a  =  (p  -  p0>  tp0 
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Other  effects,  such  as  the  Kerr  effect,  may  add  to  An.  They  will 
not  have  the  same  distribition  as  acoustic  compression,  in  general,  nor  will 
they  vary  the  same  way  with  beam  size.  However,  they  can  be  added  to  n  later. 
Let  us  ignore  them  for  the  present. 

3.3  THE  LASER  BEAM 

For  the  purpose  of  this  discussion,  the  laser  pulse  is  represented 
by  two  impulses  of  equal  energy.  Each  impulse  has  a  Gaussian  radial  intensity 
distribution.  The  radius  is  measured  to  the  point  where  the  intensity  drops 
to  1/e  of  the  peak  intensity.  The  beam  is  circularly  symmetrical  and  gently 
focused.  At  the  focus  the  intensity  distribution  is 

l(r,  t)  =  (W/ir  w1 2 3 4 5)  exp  (-r2/w2)  6  (t  ±  At/2) 

where  W  is  the  energy  of  one  impulse.  The  two  impulses  are  separated  by  a 
time,  At,  equal  to  2/3  the  half-height  duration  time  of  a  laser  pulse.  Thus 
if  the  physical  laser  pulse  has  a  Gaussian  time  distribution,  the  two  impulses 
will  occur  at  the  centroids  of  the  two  halves  of  the  physical  pulse. 

The  temporal  intensity  distribution  is  shown  in  Figure  3,  and  the 
spatial  Intensity  distribution  is  sketched  in  figure  4,  The  resulting  force 
distribution  is  given  in  Figure  5. 

3.4  THE  ACOUSTIC  WAVE 

The  electrostrictive  force  drives  n  radially  propagating  acoustic 
wave.  Usually  when  acoustic  beam  trapping  occurs  the  boundaries  of  the  mate¬ 
rial  are  so  far  from  the  beam  center  that  the  trapping  event  occurs  sooner 
than  sound  can  be  reflected  from  the  boundaries.  Hence,  the  acoustic  wave 
equation  applies  with  the  following  conditions: 

1)  Circular  symmetry 

2)  Infinite  homogeneous  isotropic  medium 

3)  Solution  at  the  beam  center  is  well  behaved 

4)  Solution  at  infinite  radius  is  zero 

5)  Solution  is  not  a  function  of  z  (paraxial  beam  case) 
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I 
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Figure  3.  Laser  Pulse  Temporal  Intensity  Distribution, 
and  the  Two-Impulse  Approximation. 
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Figure  4. 

Gaussian  Spatial  Intensity 
Distribution 


Figure  5. 

Initial  Force 
Distribution 


Figure  6. 

Initial  Rate  of  Change 
of  Compression  Distribution 


Figure  7, 

Hankel  Transform  of 
Initial  Rate  of  Change 
of  Compression 


Figure  8. 

Acoustic  Compression  at 
T  =  0.596.  Proportional 
to  Refractive  Index  Change 
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The  acoustic  wave  equation  can  be  solved  tor  the  case  of  an  impulse 
driving  force.  Since  the  equation  Is  linear,  It  Is  convenient  to  use  dimen- 
sionless  variables  and  a  normalized  impulse. 

Let  the  dimensionless  radial  coordinate  be  x,  the  radius  measured 
in  units  of  the  characteristic  beam  radius  w; 


x  =  r/w 

There  is  a  characteristic  acoustic  response  time,  T,  equal  to  the  time  re 
quired  for  sound  to  traverse  the  beam  radius; 


T  =  w/v 


where  v  is  the  speed  of  sound  for  a  two-dimensional  compression  wave, 
v  =  (Y/2 pc)1/2  (1  *  O  1/4 

Here  Y  is  Young's  modulus  and  €  is  Poisson’s  ratio  for  solid  media. 
Let  the  dimensionless  time  variable  be  T,  time  measured  in  units  of  T. 


(1) 


T  =  t/T  -  (v/w)  t 

This  choice  of  units  make's  the  velncUv  of  sound  equal  unity  in 
the  wave  equation.  Also,  the  acoustic  velocity  £  is  normalized  by  the  velocity 
of  sound,  and  the  acoustic  displacement  u  is  normalized  by  w.  The  acoustic 
compression  a  is  the  negative  divergence  of  the  displacement: 

a  =  -  V  *  U 

w'  -e  U  =  u/w.  Differentiating  with  respect  to  r,  we  have 
ba/bT  =  -  7  •  dU/dT 

where  dU/oT  =  u/v 
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The  initial  rate  of  change  of  compression  distribution  appears  in 

Figure  6. 

When  the  electrostrictive  force  acts  impulsively  on  a  medium  ini¬ 
tially  at  rest,  there  is  no  immediate  displacement  or  compression.  The  ini¬ 
tial  velocity  and  initial  rate  of  change  of  compression  may  be  deduced  from 
conservation  of  momentum, 

I7dt  "  <’o“ 

The  compression  (7  is  a  function  of  radius  and  time  only.  To  sum 
up,  the  problem  is  to  solve  the  acoustic  wave  equation: 


2 

V  CT 


n  ct/8t 


2 


subject  to 

1)  CT  =  ct(x,t) 

2)  ct(x,0)  =  0 

2  2 

3)  9ct(x, 0)At  =  A(1  -  x  )  exp  (-x  ) 


where 

A  =  (2W/3imocvpow^)  (n^  +  2)  (n^  -  1) 

4)  ct  (0,r)  4  oc 

5)  <7<bc  ,  r)  =  0 

Th"  equation  and  all  of  the  conditions  except  (3)  are  satisfied 
by  ct  (x,  t)  =  (B/A)Jo(y x)  sin  (yr)  where  B  and  y  are  arbitrary.  A  linear 
superposition  of  solutions  will  satisfy  condition  (3)  as  well.  The  correct 
cbmbination  of  solutions  g(y)  is  given  by  the  Hankel  transform  (not  the 
Fourier  transform  because  of  circular  symmetry).  The  weight  function  if  x 
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and  the  eigenfunctions  are  JQ(yx)  corresponding  to  1  and  sin  (xy)  or  cos  (yx) 
for  a  Fourier  transform.  Thus, 

m 

g(y)  -  )  i  [a(i  -  £2)  exP  (-£2)  1  J0(^d? 

=  *,y2/8)  exp  (-y2/4) 

This  distribution  is  plotted  in  Figure  1.  Now  superposition  is 
applied  to  obtain  the  solution. 

a(x,  r)  ---  A  f  (y2/8)  exp  (-y2/4)  J  (yx)  sin  (yT)  dy 
°o 

This  difficult  integral  has  been  evaluated  numerically  as  described  below. 

For  the  important  on-axis  case,  where  x  =  0,  the  integration  may  be  performed 
analytically,  since  Jq(0)  =  1« 

a(0,r)  =  A  |  (y2/8)  exp  (-y2/4)  sin  (y r)  dy 

Jo 

=  A  r/2  +  (A/2)  (1  -  2r2)  exp  (~r2)  [  exp  (£2  )  d£ 

o 

=  A  (r-  4r3/3  +  -  -  -) 

This  function  is  plotted  in  Figure  9. 

3.5  THE  NUMERICAL  INTEGRATION 

The  integral  for  q(x,  r)  was  computed  in  our  Scientific  Computer 

Facility.  We  programmed  the  Scientific  f  Systems  9300  Computer  with  65 

Fortran  IV  statements.  Simpson's  method  w«*s  used.  The  Bessel  function  was 

"8 

generated  by  a  polynomial  approximation  accurate  to  5  x  10  absolute  error 
from  the  Handbook  of  Mathematical  Functions  (Ref.  4).  One  hundred  and  one 
values  of  y  were  used  each  time  the  integral  was  evaluated.  The  integral  was 
computed  at.  101  values  of  x  for  each  of  124  values  of  j.  Time  was  saved  by 
storing  parts  of  the  kernel  which  do  not  change,  and  by  buffering  the  printer 
Running  time  was  about  40  minutes,  or  about  1  millisecond  for  each  evaluation 
of  the  Bessel  function. 
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The  computed  points  were  saved  on  magnetic  tape.  We  displayed  the 

points  with  the  on-line  oscilloscope  associated  with  our  SDS  930  Computer. 

The  Fortran  II  display  program  allowed  us  to  show  stationary  frames  continu¬ 
ously,  or  to  show  frames  sequentially  at  any  comfortable  viewing  ratt?.  The 
display  was  useful  in  studying  the  wave  motion  and  ir.  checking  the  computation. 

Refer  to  Figures  10  and  11  for  plots  of  wave  amplitude  a/A  versus 
x  at  various  values  of  T.  The  compression  at  maximum  on-axis  amplitude  is 
plotted  in  Figure  8,  for  comparison  with  the  other  functions  in  the  solution. 

3.6  THE  TRAPPING  CONDITION 

An  unexpected  result  of  the  computed  solution  of  the  acoustic  wave 

is  the  fact  that  the  compression  is  greatest  on  axis,  for  0  <  r  ^  0.85.  Thus, 

it  is,  from  the  start,  a  focusing  distribution.  There  is  no  latency  period  of 
zero  or  ne  ative  focusing  before  some  positive  focusing  begins.  The  on-axis 
solution  shown  in  Figure  9  shows  that  the  convergence  (reciprocal  focal  length) 
varies  linearly  with  7-,  for  j  «  1. 

The  two-impulse  model  is  valid  for  T  <  0.85.  The  first  impulse 
starts  an  acoustic  wave,  and  the  acoustic  wave  builds  up  a  distributed  lens 
r, f  increasing  convergence.  Diffraction  of  the  second  impulse  will  be  defeated 
if,  at  the  time  the  second  impulse  occurs, 

£n/n  >  (1/2)  (\/2imow) 

This  on-axis  An  will  be  the  sum  of  An  caused  by  the  acoustic  wave  and  Kerr 
effect.  The  acoustic  An  is  a  function  of  beam  energy  and  time  duration,  while 
Kerr  effect  depends  on  intensity.  Thus  the  two  effects  do  not  vary  the  same 
way  with  beam  size. 

3.7  THE  TRAPPING  THRESHOLD 

For  optical  glass  and  other  materials  in  which  Kerr  effect  is  weak, 
we  may  calculate  the  approximate  energy  threshold  for  trapping  by  equating  the 
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required  axial  n  with  the  \n  due  to  the  acoustic  wave.  The  total  beam 
energy  W^,  2W.  Hence 

9e\^  p  _ no  Pq _ 1  vw 

WTRAP  TT"  (n2  +  2  2  2  _  J  0/ A 

'  o  o 

For  r  '<  1,  a/ A  “*■  T  =  t  v/w 

Thu  a ,  for  pulses  that  are  short  compared  with  the  beam  radius  transit  time, 

'Jl  ,z  r _ "-I  no _ 1  H — 

™ap  4“  („2  *  2) 2  <n2  -  l)2  Lt 

o  o 

we  may  define  the  power  for  the  impulse  pair  as 
P  -  WTOTAl/“- 

Then  for  shore  duration  pulses  (Ar^l)  we  have 

P  .  [  __ - - 

TRAP  4*  L  (n‘  ,  2) 

o 

For  T  -  2'l/\  w2/At2  -  v2/2 

Thus,  when  the  pulse  duration  is  matched  to  the  beam  transit  time,  = 

P  a  Dower  level  characteristic  only  or  the  material  (and  the  wavelength). 

CRITICAL  •*  K 

The  quantity, 

P  /T 

CRITICAL' 

is  the  maximum  average  rate  at  which  power  may  be  increased  in  a  material  over 
a  time,  T,  without  acoustic  trap  formation. 


o  Po 


2  2  .  2 
(no  -  U 


1  ~2 


At 
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3.8  NUMERICAL  CALCULATIONS 

The  acoustic  trapping  threshold  curve  for  a  material  can  be  ob¬ 
tained  with  simple  calculations.  It  is  only  necessary  to  obtain  the  sound 
velocity,  the  beam  size  that  matches  the  pulse  duration  time,  and  the  critical 
power  level.  Then  simple  graphical  methods  yield  the  threshold  curve. 

Sound  Velocity:  The  relevant  sound  velocity  is  that  for  a  radially 
propagating  compressional  wave.  For  solids,  use  equation  (1)  in  Section  3.4 
The  Young's  mo  ulus,  density,  and  Poisson's  ratio  for  glasses  are  found  in  the 
newer  Schott  glass  catalogs. 

Example:  "or  BK-7,  Y  =  8310  kp/mm2,  p  =  2.51  g/cm3,  c  ■=  0.208. 

Multiply  p  by  103  to  obtain  the  density  in  kilograms  per  cubic  meter,  and  Y 

by  9.81  x  10  to  obtain  newtons  per  square  meter  from  kiloponds  per  square 

3 

millimeter.  Hence  v  =  4.27  x  10  meters  per  second. 

Matched  Beam  Size:  Let  the  pulse  duration  time  p  be  measured  at 

the  half-peak  points.  Then  w  =  pv. 

matched 


Critical  Power  Level:  The  critical  power  level  for  the  material 


is 


Critical ^8l  C^8*>  PvV  [<"o  +  2>2<"o  '  1)2  1 

Example:  For  BK-7,  at  \  =  694.3  nm,  n^  =  1.45.  The  critical 
power  level  is  440  kilowatts. 

Plotting  the  Threshold:  Use  log-log  graph  paper  having  equal  siz: 
divisions  for  the  horizontal  and  vertical  scales.  Draw  a  horizontal  line 
corresponding  to  the  critical  power  level.  Mark  the  matched  point  on  the 
critical  power  line,  corresponding  to  the  matched  beam  size. 

For  a  quick,  conservative  overestimate  of  the  trapping  threshold 
curve,  simply  draw  a  line  with  a  slope  of  2  upward  from  the  matched  point. 

(The  angle  is  63. 5°  to  the  horizontal.) 
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For  an  accurate  graph,  plot  the  curve  in  Table  I  on  a  second 
piece  of  log" log  paper.  Lay  the  first  piece  over  the  second  so  the  axes  are 
parallel  and  the  natched  point  is  just  over  the  bottom  point  on  the  second 
graph.  Trace  the  curve. 

TABLE  I 

NUMERICAL  VALUES  USEFUL  FOR 
PLOTTING  ACOUSTIC  TRAPPING  THRESHOLD  CURVES 


Horizontal  Axis 

Vertical  Axis 

0.91 

0.965 

1.00 

1.00 

1.10 

1.05 

1.25 

1.20 

1.43 

1.43 

l.b8 

1.78 

2.00 

2.36 

2.50 

3.48 

3.33 

5.90 

5.00 

12.8 

10.0 

55.0 
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SECTION  IV 


AN  ANALYTICAL  SOLUTION  OF  THE  LASER  BEAM 
RADIUS  EQUATION  FOR  KERR  EFFECT  TRAPPING 


A  high  intensity  laser  beam  passing  through  a  material  medium  can 
focus  itself  into  a  long,  thin  filament  and  propagate  without  normal  diffrac¬ 
tion  spreading.  This  self-trapping  phenomenon  arises  when  the  medium's  re¬ 
fractive  index  n  is  higher  along  the  beam  axis  than  along  the  beam  edges. 

Such  refractive  index  distribution  acts  like  a  series  of  thin  positive  lenses, 
as  shown  in  Figure  12. 

The  high  intensity  beam  can  set  up  a  focusing  refractive  index 
distribution  by  .everal  physical  mechanisms,  such  as  electrostriction,  anom¬ 
alous  dispersion,  and  reorientation  of  molecular  dipole  moments.  The  last 
effect  is  Kerr  effect,  after  its  discoverer,  John  Kerr  (1824-1907).  In 
liquids  the  molecular  reorientation  can  occur  in  times  on  the  order  of  10 
picoseconds.  The  effect  is  thus  virtually  instantaneous  compared  with  the 
duration  of  nanosecond  laser  pulses,  but  it  is  much  too  slow  to  follow  500 
t  rahertz  light  wave  fields.  The  local  change  in  index  of  refraction  is 
proportional  to  the  local  beam  intensity.  When  the  beam  power  is  above  a 
certain  threshold  power  level,  the  laser  beam  focuses  itself  to  a  smaller 
beam  size  and  higher  intensity.  The  smaller,  higher  intensity  beam  causes 
still  stronger  self-focusing  until  the  beam  has  trapped  itself  at  a  small 
radius,  limited  by  diffraction. 

4.1  TUB  BEAM  RADIUS  EQUATION 

This  note  sets  forth  an  analysis  of  Kerr  effect  trapping.  The 
analysis  is  based  on  a  ray  tracing  equation  published  by  Tien,  Gordon,  and 
Whinnery  (Ref.  6).  Equation  (9)  in  that  article  reduces  to 

32a/?tz2  =  a  3  (\/?irno)2  +  (a/nQ)  }2n(0,  z,  t) /^r2  (2) 
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where 


a(z,t)  =  beam  radius  to  the  point  where  the  intensity  is 

1/e  of  the  peak  intensity, 

n(r,z,t)  =  local  index  of  refraction 

n  =  undisturbed  or  nominal  index  of  the  medium 
o 

\  =  vacuum  wavelength  of  the  laser  beam, 
r,2  =  cylindrical  coordinates  for  the  beam. 


For  this  equation  to  be  valid  Tien,  Gordon,  and  Whinnery  require 

that 

1)  "the  light  beam  is. .. launched  with  a  Laguerre-Gaussian 
or  Hermite-Gaussian  field  distribution" 

2)  "the  refractive  index  of  the  medium  varies  slowly  in  space 
(negligibly  in  an  optical  wavelength)" 

3)  the  variations  in  n  are  small  compared  with  nQ 

The  last  two  restrictions  are  violated  when  the  beam  is  trapped 
to  a  filament  with  approximately  one  wavelergti.  ^adius,  and  the  scattered 
light  is  spread  over  a  wide  wavelength  range.  However,  the  equation  can  be 
used  to  study  the  collapse  of  the  beam  toward  the  trapped  condition,  and  to 
study  the  phenomenon  when  the  beam  power  is  below  threshold. 


4.2  BEAM  PROPAGATION  WITHOUT  SELF-FOCUSING 


When  the  beam  intensity  is  weak,  the  last  term  in  equation  (2) 
is  negligible.  Thu  ray  paths  are  then  hyperbolas.  The  beam  radius  solution 
is 


a(z)  =  (a 


(\z/2rtn  a  ) 
o  o 


!] 


1/2 


(3) 


where  a  and  a.  are,  respectively,  the  initial  radius  and  initial  slope  of  the 
o  1 

rays  just  after  entering  the  medium  at  z  =0.  For  example,  let  X  =  lum, 

n  =  1.5,  a  =  0.1mm,  a,  =  -0.001839.  Then  the  beam  is  focused  to  a  radius  of 
o  o  1 

0.05mm  at  z  =  4.1cm.  This  is  the  case  of  a  very  gently  focused  beam.  The 


24 


PERKIN  ELMER 


Report  No.  8995 


first  frame  of  Figure  14  shows  a  plot  of  a(z)  for  this  case  .  Note  that  the 
vertical,  radius  scale  is  greatly  exaggerated. 


4.3  THE  LASER  BEAM  INTENSITY  EQUATION 

The  lowest  order  mode  for  the  laser  beam  has  a  nearly  Gaussian 
intensity  profile  characterized  by  the  radius,  a,  to  the  point  where  the 
intensity  drops  to  1/e  of  the  peak  intensity.  If  the  propagation  medium  is 
not  too  violently  inhcmogeneous  and  if  the  inhomogeneity  is  radially  sym¬ 
metrical,  the  Gaussian  profile  is  maintained  along  the  beam,  although  the 
beam  radius  varies  irec^use  of  diffraction  and  focusing  by  the  inhomogeneities. 

The  laser  power  may  also  vary  as  a  function  of  time.  For  Q-switched 
pulse  lasers  generally  the  pulse  energy  can  be  measured  and  some  rough  idea  of 
the  time  distribution  of  the  pulse  can  be  obtained.  For  purposes  of  analysis 
we  can  use  a  simple  pulse  shape  such  as  the  unit  quartic  pulse  shown  in 
Figure  13.  This  shape  has  a  continuous  derivative  and  finite  extent. 

Thus  the  laser  beam  intensity  is 


I(r ,  z,  t)  =  (W/ira^)  exp  (-r  /a  )  q  (t;p) 
where  W  is  the  pulse  energy  and  q(t;p)  is  the  unit  quartic  pulse  of  time 


constant  p: 


;;p)  *  [o 


(15/16p)  (t/p) 2(2-t/p) 2,  0  s  t  s  2p 


otherwise. 


Of  course,  a  is  a(z,t)  given  by  the  solution  of  equation  (2).  The  intensity 
distribution  is  so  normalized  that 


j>r 


r  dr  I  dt  I  (r,z,t)  =  W 


independent  of  the  value  of  z. 
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4.4  VARIATION  OF  REFRACTIVE  INDEX  WITH  INTENSITY 


The  change  in  refractive  index  can  be  expressed  in  cgs  units  as 

2  ,  XT-2 

n  =  n  +  —  \JE  +  . 

i 

where  J  is  the  high  frequency  Kerr  constant  due  to  molecular  rotation.  Thus 
the  change  in  index  is  proportional  to  the  local  intensity  for  either  a  plane 
polarized  or  circularly  polarized  beam  (although  the  constants  of  proportion¬ 
ality  differ  for  the  two  cases).  Let  the  constant  of  proportionality  be  K, 
depending  on  the  wavelength,  the  material,  and  the  beam  polarization. 

n(r,z,t)  =  nQ  +  KI  (r,z,t) 

In  equation  (2)  the  last  term  becomes 

-a"3(2  KW/im  )  q  <t;p) 
o 

Since  both  terms  on  the  right  of  equation  (2)  are  proportional  to 
a  3,  the  solution  at  any  time  will  be  a  hyperbola  whose  vrlue  and  slope  at 
0  are  aQ  and  a^,  respectively.  However,  the  focal  point  and  semiaxes  of  the 
hyperbola  may  vary  with  time.  The  equation  becomes 

2  2  -3 

d  a/3z  =  f(t)a  J 


where 

f(t)  =  (\/2jmo)Z  -  2KWq(t;p)/nnQ. 

The  (solution  is 

r  2  -2  2  T 

a(z, t)  =  [(aQ  +  aLz)  +  f(t)  aQ  z  J 

A  negative,  zero  or  complex  value  of  a  would  be  a  non-physical 

solution.  This  places  a  restriction  on  the  minimum  value  of  f(t)  which  occurs 

2 

when  t  *  p.  We  have  a  (z,p)  equal  to  a  quadratic  expression  in  z.  In  order 

o 

that  a  (z,p)  be  greater  than  zero  for  all  real  values  of  z,  the  discriminant 

2 

of  the  quadratic  expression  must  be  negative.  Then  the  only  zeros  of  a  will 
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be  complex.  This  condition  implies  that  £(t)  must  always  be  positive.  The 
minimum  value  of  f(t)  occurs  when  t-p,  i.e.,  at  the  peak  of  the  laser  pulse. 
For  a  valid  solution  we  must  have 


1  2 
£n/nc  <  2  (X/2nnQa) 

where  £n,  the  on-axis  increase  in  refractive  index,  is 

2 

£n  =  KWq(p;p) /na 


This  condition  is  equivalent  to  the  Pierce  stability  criterion  L/f  <  4.  where 
L  and  f  are  the  spacing  and  focal  length  of  lenses  in  a  sequence  (Ref.  5). 


The  threshold  value  of  KW/p  becomes  (KW/p) 

threshold 


Let  k  be  the  fraction 


2\2/rtn 


o 


k  =  (KW/p) /(KW/p) thre5hold 


Then  f(t)  becomes 

f(t)  =  (\/27tnQ)2  [l-kq(t;p)/q(p;p)j 

Refer  to  Figures  14  through  16  for  plots  of  the  beam  radius  when  k  =  0.999, 
0.95,  and  0.85. 
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SECTION  V 

BEAM  TRACING  IN  THE  ACOUSTIC  TRAP 


The  theory,  developed  from  the  two-impulse  model  in  Sections  II 
and  III,  yields  a  threshold  for  acoustic  beam  trapping  which  is  below  the  ex¬ 
perimental  damage  thresholds  for  the  glassy  materials  tested.  It  was  solved 
by  completely  analytical  means,  using  tabulated  properties  of  special  functions. 
The  computer  was  only  used  to  obtain  values  for  the  radial  sound  wave  off  the 
axis  of  the  beam-  That  computation  was  not  essential  to  the  solution,  although 
it  did  give  the  following  Important  insight:  The  refractive  index  dis:ribution 
set  up  by  the  acoustic  wave  is  initially  a  focusing  distribution,  al though  its 
strength  varies  with  time.  However,  the  theory  did  not  allow  for  the  continuous 
interaction  of  light  and  sound,  nor  did  it  show  how  the  focusing  action  develops 
with  time. 

5.1  OTHER  FEATURES  OF  THE  GLASS  DAMAGE  PHENOMENON 

A  complete  theory  should  also  explain  the  other  observed  features 
of  internal,  filamentary  glass  damage,  such  as  the  following: 

1)  The  Spectrum  of  the  Side-Scattered  Light 

A  white  light  flash  is  seen  when  the  damage  event  occurs. 

0  O 

Ruby  laser  light  may  be  shifted  from  6943A  down  to  4000A. 

Such  large  shifts  can  only  be  explained  by  a  very  strong 
dynamic  refractive  index  change. 

2)  Starting  Location  of  the  Damage  Track 

Even  if  the  laser  is  focused  at  the  entrance  face  of  the 
glass  samole,  the  damage  track  never  begins  at  the  entrance 
face.  There  is  always  a  short  interval  between  the  entrance 
face  of  the  sample  and  the  start  of  the  track. 
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3)  Exit  Face  Pitting 

Usually  the  track  ends  in  a  pit  on  the  exit  face  of  the 
glass  sample.  As  a  matter  of  fact,  the  threshold  for  the 
exit  face  pitting  appears  to  be  a  little  lower  than  the 
threshold  for  track  formation.  Nevertheless,  the  phenomena 
is  different  from  the  type  of  surface  damage  reported  by 
others. 

4)  Location  of  Damage  Stars 

The  damage  track  is  not  always  continuous;  it  may  start  and 
stop  several  times  in  the  sample.  Often  there  is  a  damage 
star  on  the  upstream  end  of  the  track.  These  damage  stars 
are  localized  regions  of  gross  fracture.  Occasionally 
they  show  discoloration,  indicating  possible  chemical  de¬ 
composition.  If  the  damage  track  extends  upstream  of  the 
damage  star  it  usually  only  does  so  for  a  very  short  distance 
compared  to  the  extension  downstream  from  the  damage  star. 

The  location  of  the  damage  star  is  intrigu.'ng.  If  the 
damage  starts  in  the  damage  star  and  then  propagates  down¬ 
stream  to  form  the  track,  why  doesn't  the  damage  star  cast 
a  downstream  shadow  and  prevent  the  beam  from  concentrating 
in  the  thin  filament? 

5)  Track  Propagation  Speed 

Since  the  damage  event  occurs  in  nanoseconds,  its  dynamics 
are  difficult  to  follow.  There  is  experimental  evidence 
that  the  event  which  forms  the  track  propagates  at  aiout 
10  times  the  speed  of  sound  in  the  glass.  The  experimental 
evidence  does  not  show  unequivocally  that  the  propagation 
is  either  upstream  or  downstream. 

Most  of  these  features  are  explained  qualitatively  in  the  theory  of 
electrostrictively  driven  acoustic  beam  trapping. 
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5.2  DEFINITION  OF  THE  PROBLEM 

Our  purpose  is  to  understand  how  the  acoustic  beam  trap  forms,  not 
how  the  damage  occurs.  Therefore  we  only  need  consider  small  acoustic  pressure 
changes  in  the  material.  We  will  see  how  these  lead  to  a  beam  instability  when 
the  oower  exceeds  a  certain  threshold.  The  nature  of  the  instability  is  such 
that  the  beam  rapidly  focuses  itself  into  a  thin  filament.  The  damage  occurs 
when  the  beam  is  in  this  trapped  condition.  The  instability  itself  can  be 
studied  by  means  of  a  geometric  r~v  tracing  equation  which  includes  first  order 
diffraction  effects.  This  beam  tracing  equation  is  equation  (2)  of  Section  IV. 


The  sound  wave  satisfies  essentially  the  same  conditions  as  those 
given  in  Section  III,  except  that  it  is  driven  continuously  by  the  light  wave, 
rather  than  running  inertially  after  an  impulse.  For  a  gently  focused  beam, 
axial  components  of  the  sound  wave  field  are  negligible  in  comparison  with 
the  radial  components. 

The  desired  solution  will  show  a  graph  of  the  beam  radius  plotted 
versus  axial  length  for  each  of  many  small  time  intervals  during  the  laser 
pulse.  The  time  distribution  of  the  laser  pulse  can  be  assumed  to  be  the 
quartic  pulse  distribution  described  in  Section  IV,  since  the  actual  details 
of  the  shape  of  the  laser  pulse  do  not  have  too  much  bearing  on  the  threshold 
for  damage  and  on  the  tire  development  of  the  trap. 


5.3  THE  COUPLED  EQUATIONS 


The  sound  wave  equation  is 


2  2 

iW  + 


r_1ao/dr 


2  .  2 
b  o/?2 


-22.2 

v  5 

-  (l/6vn^o)  (nQ  +  2)  (n^  -  1)  (W/jt)q(t  ;p) 

.  -4,,  2.2.  .  2,2. 

•  4  a  (1-r  /a  )  exp  (-r  /a  ) 


The  beam  tracing  equation  is 

d2a/gz2  =  a‘3(K/2rtno)2  +  (a/6n2)(n2  +  2)  (n2  -  1)  d2(j(o,  z,  t)  /*r2 
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The  general  solution  of  these  coupled,  nonlinear,  partial  dif¬ 
ferential  equations  has  been  left  as  an  exercise  for  the  computer.  One 
simplification  can  be  made  immediately.  The  beam  tracing  equation  involves 
only  the  second  derivative  with  respect  to  r  of  the  on-axis  sound  wave  solution. 
Since  che  sound  wave  equation  is  linear  we  Tiay  take  its  Hankel  transform  ana¬ 
lytically.  When  z  is  held  fixed,  the  Hankel  transform  of  (j  satisfies  an 
ordinary  differential  equation  in  t.  This  ordinary  equation  may  be  solved  using 
initial  values  of  the  Hankel  transform,  and  of  its  derivative  with  respect  to 
time,  and  assuming  a  step-plus-ramp  driving  function.  The  solution  is  then 
found  by  taking  the  inverse  Hankel  transform.  Since  the  beam  tracing  equation 
involves  only  the  second  derivative  of  the  sound  wave  amplitude  at  r  =  0,  and 
the  inverse  Hankel  transform  involves  r  only  as  part  of  the  argument  of  the 
zeroth-order  Bessel  function,  we  may  perform  the  differentiation  analytically 
under  the  integral  sign.  Thus  for  each  time  step  and  each  value  of  z,  the 
sound  wave  equation  can  be  solved  with  only  one  numerical  integration.  The 
use  of  a  step-plus-ramp  driving  function  allows  us  to  take  a  much  coarser 
time  step  than  we  would  require  if  we  used  an  impulsively  driven  equation. 

5.4  THE  STEP-BY-STEP  COMPUTER  SOLUTION 

The  problem  is  divided  up  into  sections  and  slices  as  shown  in 
Figure  17.  In  its  present  state  of  development  the  computer  program  can  hanu^e 
101  sections.  The  zeroth  section  is  the  entrance  face  of  the  solid  ov  liquid 
medium. .  The  initial  beam  radius  and  entrance  angle  are  fixed  on  the  zeroth 
section.  Since  the  solution  is  radially  symmetrical,  it  will  be  the  same  in 
each  radial  slice. 

The  computer  program  must  store  values  of  a  for  each  of  the  101 
sections.  Also  for  each  section,  it  must  store  40  values  of  the  Hankel  trans¬ 
form  of  the  sound  wave,  and  40  values  of  the  derivative  of  the  Hankel  transform 
of  the  sound  wave  with  respect  to  time.  An  array  of  101  points  is  also  set 
aside  for  the  second  derivative  of  the  sound  wave  with  respect  to  r.  No 
Bessel  function  storage  or  Bessel  function  subprogram  is  required.  Initially 
the  Hankel  transform  of  the  sound  wave  and  its  derivative  with  respect  to  time 
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Figure  17  •  Geometry  of  the  Beam  Trapping  Analysi 
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are  set  to  zero  along  with  the  second  derivative  of  the  sound  wave  with  re¬ 
spect  to  r.  The  solution  then  proceeds  in  steps: 

1)  Solve  the  beam  tracing  equation  using  the  Runge-Kutta  method. 
The  axial  step  size  in  the  solution  is  1/5  of  the  interval 
between  sections.  Values  of  the  second  derivative  of  the 
sound  wave  with  respect  to  r  are  interpolated  by  Lagrange 
cubic  fitting. 

2)  Store  101  values  of  the  beam  radius  a.  Also  save  these 
values  on  tape  for  later  plotting. 

3)  If  the  time  has  not  yet  reached  2. Op,  the  end  of  the 
laser  pulse,  advance  the  time  by  a  snail  time  step  and 
obtain  the  new  value  of  the  beam  power  at  that  time. 

4)  For  each  section  of  the  beam  compute  the  second 
derivative  of  the  sound  wave  with  respect  to  r  by 
using  Simpson  intergration.  At  the  same  time  revise 
the  value  of  the  Hankol  transform  of  the  sound  wave 
and  its  time  derivative  for  the  next  time  step. 

5)  Go  to  1. 

The  program  parameters  are  the  wavelength,  the  nominal  refractive 
index,  the  initial  beam  radius  and  the  initial  sW  r r  .  beam  radius  at 
the  entrance  face,  and  the  fraction  of  pulse  power  divided  by  theoretical 
threshold. 

5.5  TYPICAL  RESULTS  FOR  ACOUSTIC  BEAM  TRAPPING 

Figure  18  depicts  acoustic  beam  trapping  when  the  pulse  power  is 
3007,  of  the  threshold.  There  are  11  frames  showing  the  beam  radius  plotted 
versus  axial  length  ror  11  different  times  during  the  pulse.  The  first  frame 
shows  the  path  of  the  beam  when  the  illumination  has  just  begun.  The  sixth 
frame  gives  the  trace  when  the  pulse  power  is  highest.  The  11th  frame  shows 
the  path  taken  by  the  light  in  the  trailing  end  of  the  pulse. 
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Figure  18 .  Acoustic  Beam  Trapping  at  300  Percent  Threshold 
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The  action  may  be  interpreted  by  considering  that  each  section  of 
the  beam  forms  an  acoustic  lens.  Initially  the  lenses  have  zero  power,  but 
some  of  them  build  up  focusing  power  faster  than  others.  The  strongest  lens 
is  formed  at  the  focus  of  the  beam.  As  ii_s  focusing  power  builds  up  in  time 
it  causes  the  beam  to  be  refocused  at  some  distance  downstream.  As  the 
strength  of  the  acoustic  lens  continues  to  build  up,  its  focal  length  is 
shortened  so  the  second  focus  moves  upstream.  The  second  focus  likewise  forms 
a  third  focus  and  so  on.  A  little  after  the  laser  pulse  is  half  over,  three 
definite  foci  have  appeared  in  the  12-l/2cm  region  plotted.  At  the  same  time 
that  these  acoustic  lenses  are  becoming  stronger  because  of  the  inertial 
properties  of  the  material,  the  laser  pulse  power  is  decreasing.  The  net 
result  is  regions  of  sharper  and  sharper  focus  interspersed  by  regions  in  which 
the  beam  radius  is  increasing.  Eventually  the  beam  will  be  able  to  escape 
again  to  large  radii.. 

Note  in  the  last  frame  that  the  sharpest  focus  occurs  the  farthest 
downstream.  This  is  to  be  expected  since  the  focus  there  is  the  result  of 
two  strong  lenses  upstream.  Another  important  feature  is  the  motion  of  the 
first  focus.  Notice  how  it  moves  upstream  in  the  beam.  This  is  because  the 
acoustic  lenses  in  the  first  few  sections  of  the  sample  are  also  developing 
in  strength  with  time. 

Several  features  of  the  damage  phenomenon  are  immediately  ex¬ 
plained.  The  fact  that  white  light  is  scattered  largely  from  the  region  of 
the  damage  stars  is  easily  seen  to  follow  from  the  fact  that  the  highest 
focusing  and  therefore  the  highest  index  of  refraction  occur  at  places  which 
become  damage  stars.  Also  the  track  propagates  backward,  becoming  more  and 
more  sharply  defined  until  it  ends  in  a  damaged  star.  If  there  are  more  than 
one,  the  damage  stars  occur  in  backward  sequence.  The  farthest  downstream 
occurr  first,  then  the  next  one  upstream,  then  finally' the  damage  star  closest 
to  the  upstream  end  of  the  beam.  Thus,  the  damage  event  does  not  cast  a  shadow 
that  gets  in  its  own  way. 
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The  reason  why  it  is  difficult  to  measure  the  spee:’  of  propagation 
of  the  damage  event  is  now  apparent.  The  intensity  varies  widely  and  more  than 
one  focus  region  may  move  past  the  velocity-sensing  optics. 

The  pit  on  the  exit  face  of  the  glass  sample  is  formed  by  axial 
sound  wave  components  which  do  become  large  when  the  beam  is  sharply  focused. 
There  can  never  be  a  damage  pit  on  the  entrance  face.  The  entrance  conditions 
of  the  beam  are  fixed;  therefore,  the  bean  cannot  focus  itself  sharply  at  the 
entrance  face  of  the  glass. 

All  of  the  qualitative  features  of  the  phenomenon  are  a  sharp  con¬ 
trast  to  the  phenomenon  of  Kerr  effect  trapping,  in  which  there  is  one  focus 
region  which  moves  downstream. 
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SECTION  VI 

SUMMARY  OF  RESULTS  AND  FUTURE  PLANS 

During  the  first  half  of  the  contract,  we  have  accomplished  the 

following: 

1)  Kerr  Effect  Trapping  Analyzed;  An  analytical  solution  was 

** 

obtained  for  Kerr  effect  trapping.  Computer  plots  were 
made  showing  the  time  development  of  this  trapping  phenomenon. 

A  threshold  for  the  effect  was  obtained. 

i 

| 

2)  Two-Impulse  Model  Improved;  Some  corrections  and  improve¬ 
ments  were  made  in  the  earlier  two-impulse  model  for  the 
acoustic  beam  trapping  phenomenon. 

3)  Acoustic  Trapping  Equations  Derived:  Equations  were  derived 
for  acoustic  beam  trapping  when  the  light  pulse  is  treated 
as  having  a  finite  time  distribution. 

4)  Driven  Sound  Wave  Programmed:  Programs  were  written  and 
tested  to  show  the  development  of  the  sound  wave  with  a 
time  distributed  forcing  function. 

5)  Beam-Tracing  Program  Completed:  Programs  were  written  to 
compute  and  display  the  acoustic  beam  trapping  phenomenon 
in  which  the  sound-light  interaction  is  carried  out. 

6)  Combination  Program  Started;  A  computer  program  to  study 
trapping  phenomena  when  both  Kerr  effect  and  acoustic  trap¬ 
ping  are  operative  has  been  written  and  partly  tested. 


41 


i»i  «,*«. .  11  n  •  .tffliir  ,»;::"'!uiniii:..,;|4',.jiiww  jam-  ■ , HmtntiBimBrytmi i ^mmiH4ims«94*^^raa39i^MiJltKK(l!yBai)li^!C3t!^!6^l!W,tK)#S8iattl,Bnjtt!to4(,tittHS£3tttaiBaBfl8flSfS!^l9WIJ! I ^ liii 


PERKIN-ELMER 


Report  No.  3995 


[ 

l 

l 

[ 

[ 

[ 

l 

c 

[ 

r 

E 

r 

L 

L 

i 

L 

i 

L 

L 


During  tht  second  half  of  the  contract,  we  expect  to  complete 
the  development  of  the  combination  Kerr  effect  and  acoustic  trapping  program. 

Ue  will  prepare  an  article  for  publication  to  set  forth  the  theory  and  dis¬ 
play  several  cases  of  interest.  If  rime  permits,  we  will  also  prepare  a 
companion  article  which  sets  forth  the  experimental  results  from  our  old  glass 
damage  study  contract  in  more  detail.  Another  possibility  will  be  to  study 
the  reduced  crapping  thresholds  expected  if  the  material  already  has  a  focusing 
index  distribution  in  it  due  to  thermal  strains. 

A  second  quarterly  letter  report  and  a  final  report  will  be 

submitted. 
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selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  mill- 
tary  project  code  name,  geographic  location,  may  be  used  as 
key  words  bat  will  be  followed  by  an  indication  of  technical 
context.  The  assignment  of  links,  rules,  ond  wc»tgh.s  is 
optional. 


45 


